On the action potential as a propagating density 
pulse and the role of anesthetics 
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The Hodgkin-Huxley model of nerve pulse propagation relies on ion currents through specific resistors called ion channels. 
We discuss a number of classical thermodynamic findings on nerves that are not contained in this classical theory. Particularly 
striking is the finding of reversible heat changes, thickness and phase changes of the membrane during the action potential. Data 
on various nerves rather suggest that a reversible density pulse accompanies the action potential of nerves. Here, we attempted 
to explain these phenomena by propagating solitons that depend on the presence of cooperative phase transitions in the nerve 
membrane. These transitions are, however, strongly influenced by the presence of anesthetics. Therefore, the thermodynamic 
theory of nerve pulses suggests a explanation for the famous Meyer-Overton rule that states that the critical anesthetic dose is 
linearly related to the solubility of the drug in the membranes. 



1 Introduction 

The description of electrical phenomena in nerves is among the 
first biological problems studied in physics. Galvani no- 
ticed that the legs of dissected frogs made active movements 
when their nerves were connected to a battery. He called this 
phenomenon "animal electricity". After learning about these 
experiments, Volta |2] stated that nerve pulses are electrical 
conduction phenomena. Helmholtz performed the first mea- 
surements of the propagation velocity of nerves and found a 
value of about 30 m/s in the nerves from frog muscle. In the 
second half of the 19th century Ostwald and others devel- 
oped the theory of osmosis and electrochemistry, and attempts 
were made to relate the flux of ions through the nerve mem- 
branes to the propagating action potential 0. This finally re- 
sulted in the model by Hodgkin and Huxley from 1952 
that is the presently accepted model for the nerve pulse. This 
model relies on ionic currents through ion-selective objects (ion 
channel proteins) and the membrane capacitor. In the con- 
text of their model, the conductance of these objects displays 
rather complex voltage and time dependences that enter the 
differential equation via a set of empirical parameters. Those 
parameters are taken from experiment but do not yet have a 
satisfying theoretical justification. Even though Hodgkin and 
Huxley did not originally specify the ion-conducting ob- 
jects, it was clear from the line of argument that these objects 
were expected to be specific proteins called ion-channels. In 
1976, Neher and Sakmann using the patch clamp technique de- 
scribed such channels microscopically 0. Nowadays, many 
investigators all over the world investigate the properties of ion 
channels. In 1998, MacKinnon and collaborators crystallized 



the potassium channel and suggested a pathway for the potas- 
sium through a pore within the protein |8j. Thus, the Hodgkin- 
Huxley model seemingly finds support in independent experi- 
ments. The model by Hodgkin and Huxley is a purely electrical 
description based on conductors (ion channels and the cytosol 
of the nerve axon) and on a capacitor, which is the lipid mem- 
brane. It does not contain any thermodynamical variable ex- 
cept the membrane potential. Entropy, temperature, pressure 
and volume do not play a role. There is, however, strong evi- 
dence that phenomena during the action potential are not purely 
electrical. It has been observed by a number of investigators 
that the dimensions of the nerve change in phase with voltage 
changes and that the nerve exerts a force normal to the mem- 
brane surface l9l 1101 ITT1 1121 IT3l . Further, during the action 
potential lipid membrane markers change their fluorescence in- 
tensity and their anisotropy I14II15I . Most striking, however, 
is the finding that there are reversible changes in temperature 
and heat during the action potential (l6|[T7l[T8][T9l|20|. While 
the Hodgkin-Huxley model |<6] contains resistors that should 
generate heat during the flow of ions, the reversible release 
and re-absorption of heat does not find a satisfactory explana- 
tion within this model 1211 . Recently, Heimburg and Jackson 
1 22 23 1 proposed that the action potential is rather a propagat- 
ing density pulse (soliton), and therefore an electromechanical 
rather than a purely electrical phenomenon. This corresponds 
to a localized piezoelectric sound pulse within the nerve mem- 
brane. Such a model is able to explain most of the thermody- 
namical findings on nerves and results in the correct propaga- 
tion velocity of about 100 m/s for a myelinated nerve. Interest- 
ingly, Hodgkin and Huxley themselves proposed the possibil- 
ity that the nerve pulse is a propagating mechanical wave 1241 . 
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Figure 1: Left: Action potential adapted from the original paper of Hodgkin and Huxley $6\/. Right 
top: Electrical currents in the Hodgkin- Huxley model through ion channels. Right bottom: Equiva- 
lent circuit picture replacing ion channels by resistors and the membrane as a capacitor. 



Anesthesia is a phenomenon that seems to be closely related to 
the action of nerves. Since the standard model of nerve action 
is based on the action of ion channels, most research has been 
dedicated to investigating the influence of anesthetics on such 
proteins. However, an old finding by Meyer 1251 and Overton 
1261 1271 states that the action of anesthetics is linearly related 
to their solubility in membranes. This includes the noble gas 
Xenon. Although some ion channels are influenced by some 
anesthetics, there is no quantitative correlation with the well- 
documented Meyer-Overton rule 1 28 1. 

In this paper we briefly discuss some of the historical find- 
ings on nerves, including the Hodgkin-Huxley model and ther- 
modynamic data on nerves. It is shown that the Hodgkin-Huxley 
theory does not describe the thermodynamics of the nerve pulse 
correctly. Instead, the propagation of a density pulse is shown 
to explain in a quantitative manner many features of the nerve 
pulse, including density, fluorescence anisotropy and heat change 
Finally, we show that such a description leads to a satisfactory 
quantitative explanation of general anesthesia. 

2 The Hodgkin-Huxley model 

In the Hodgkin-Huxley model |6) the propagation of a volt- 
age pulse is the consequence of ion currents through the mem- 
brane and along the nerve axon. The electrochemical potential 
(Nernst potential) across the nerve membrane balances the ion 
concentration differences on both sides of the nerve axon. The 
transient opening of voltage-dependent ion channels leads to a 
related transient voltage change that can propagate. Most of the 
data on which the Hodgkin-Huxley model is based originate 
from voltage-clamp experiments on giant squid axons where 
the trans-membrane voltage is kept constant along the whole 
length of the axon. 

The relation for the ion current through the membrane un- 
der voltage clamp conditions is based on an equivalent circuit 



picture that is schematically shown in Fig. 1. Describing ion 
channels by resistors and the membrane as a capacitor, one ob- 
tains 

I m = C m ^-+g K (U-E K )+g Na {U-E Na )+g L (U-E L ) 
at 

(1) 

where I m is the current through the membrane, and C m is 
the capacitance of the membrane (typically on the order of 
ImF/cm 2 ). The Ek , E^ a and El are resting potentials that 
depend on ion concentrations. The qk and g^ a are the con- 
ductances of K-channels and Na-channels, and g^ describes the 
leakage currents. The conductances are not constants but rather 
complicated functions of time and voltage, gx = gi<(V, t) and 
9Na = gNa(V, t), that have been empirically fitted by Hodgkin 
and Huxley |:6| using many ad hoc parameters. Therefore, the 
seemingly simple eq. Q is in fact very complicated, and all the 
'mysteries of the observed phenomena are hidden in the func- 
tional dependences of the conductances on time and voltage. 
The trans-membrane current in eq. (1) is given as the sum of a 
capacitive current and an Ohmic current. The capacitive current 
is given by 

Ic = ^(C m -U) = C m ^ + U^ (2) 
at at at 

A closer look at the right hand side of eq. Q indicates that 
the capacitive current used by Hodgkin and Huxley consists 
only of the C m ■ dU /dt term and that the capacitance C m was 
assumed to be constant. Therefore the U ■ dC m /dt term has 
been neglected. This is probably not correct since we will show 
in the next section that the thickness of nerves changes during 
the pulse. Note in particular that the function dCm/dt carries 
the same units as the conductances, gi. For this reason it may 
not always be trivial to distinguish currents through resistors 
and capacitive currents in an experiment during a propagating 
pulse 1291 1301 . To arrive at a wave equation for the nerve axon, 
Hodgkin and Huxley assumed that the total current is the sum 



2 



Km Na 




Q-=HJ D-=H] \Y=ri} — □-=-□ 



pulse 



Figure 2: Equivalent circuit picture of a propagating voltage pulse. Currents flow along the nerve 
axon and across the axonal membrane through resistors and should produce net heat dissipation. 



of the trans-membrane current and the current along the axon. 
A further ad-hoc assumption is that a propagating solution ex- 
ists that fulfills a wave equation. Hodgkin and Huxley6 arrived 
at the following differential equation for the propagating nerve 
pulse: 



d 2 U 



2Ri dx 2 



E k ) + 9Na(U - E Na ) (3) 



where a is the radius of the axon and Ri is the resistance of 
the cytosol within the nerve. This equation introduces a de- 
pendence of the pulse propagation on the nerve radius. The 
elements of the propagating pulse are summarized in Fig. 2 
that shows the equivalent circuits as an in-line arrangement of 
many local equivalent circuits as shown in Fig. 1. Due to the 
voltage and time dependence of the conductances in eq. this 
differential equation can only be solved numerically. Hodgkin 
and Huxley found a convincing agreement between the calcu- 
lated and the observed pulse shape for the squid axon that only 
contains K- and Na-currents. 

One immediate implication of the Hodgkin-Huxley model 
is that ion currents through the nerves should produce heat. 
Electrical currents through resistors generate heat, independent 
of the direction of the ion flux. The heat production in such an 
experiment therefore should always be positive if the Hodgkin- 
Huxley model is taken seriously and the analogy of ion currents 
through protein pores and Ohmic currents is assumed to be cor- 
rect. The heat dissipation should be related to the power of a 
circuit through the resistor, i.e. dQ/dt = P = U-I = g~ x I 2 > 
for each of the conducting objects in all phases of the action 
potential. In the next section we will show that this is not in 
agreement with the experiment. 



3 Thermodynamics of nerve pulses 

The Hodgkin-Huxley model 1 6 1 is a purely electrical theory. It 
is based on equivalent circuits and makes use of capacitance, 
resistors and ionic currents. It is not a thermodynamic theory. 
It does not explicitly contain temperature and heat or other ther- 
modynamic variables such as pressure, volume and the chem- 



ical potentials of molecules dissolved in the membrane (e.g. 
anesthetics). However, there are many reports in the literature 
indicating that, in addition to the electrical response of nerves, 
other variables also change, for example the thickness, the en- 
thalpy and heat content of the nerve. In the following we briefly 
discuss some of these data. 

3.1 Thickness and forces 

I. Tasaki and collaborators have published several studies on the 
mechanical and thermodynamic properties of various nerves 

|9][T0|[ni[I2[I3m[m|20||3ll- For a11 nerves investigated, 
they found that the action potential (i.e. the voltage pulse) is ac- 
companied by changes in the dimensions of the nerve. In Fig. 3 
(right) it is shown that the voltage pulse of a squid axon is ex- 
actly proportional to the change of its thickness l9l l!0l . In the 
example this thickness change is about 1 nm. Further, the same 
authors showed that during this pulse a considerable force acts 
on a piston that was brought into contact with the nerve surface 
(Fig. 3, left). The force on that piston (0.01 cm 2 cross section) 
was shown to be about 2 nN at the voltage peak maximum. 

3.2 Fluorescence changes, optical changes and 
alterations in lipid state 

During the action potential not only thickness and pressure on 
a piston change but also the state of the membrane as measured 
by the fluorescence changes of lipid dyes. Tasaki and cowork- 
ers 1 1 41 1131 found that in various nerves under the influence of 
the action potential the fluorescence intensity change is propor- 
tional to the voltage pulse (see Figure 4). 

In the same paper they showed that the fluorescence ani- 
sotropy of these markers also changes (data not shown). The 
fluorescence anisotropy is a measure of the rotational mobil- 
ity of the fluorescence markers. A lower anisotropy indicates 
faster movement, whereas a high anisotropy indicates slower 
movement. Since the fluorescence anisotropy changed dur- 
ing the voltage pulse, Tasaki and collaborators 1141 concluded 
that the viscosity of the membrane changes during the nerve 
pulse. Note that they published this paper prior to the fluid 
mosaic model by Singer and Nicholson |32| from 1972 that 
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Figure 3: Mechanical changes during the action potential. Left: Force on a piston during the action 
potential in a squid axon. The solid line represents the voltage changes, the dotted curve the force. 
Right: During the nerve pulse in a squid axon the thickness of the nerve changes proportional to the 
voltage. Data on squid axons adapted from ref. $9$. 
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Figure 4: Voltage changes (top traces) and fluorescence changes (bottom traces) for 4 different 
fluorescence markers and nerve preparations. They are exactly in phase. 1. Squid giant axon and 
8-anilinonaphtalene-l -sulfonate (ANS). 2. Crab leg nerve with fluorescein isothiocyanate (FIT). 3. 
Squid axon with FIT. 4. Crab leg nerve with lysergic acid diethylamide (LSD). From ref. M4V with 
permission. The data were taken as a proof for changes of the viscosity within the membrane during 
the action potential. 



established the present view of the biological membrane. The 
concept of phase transitions in lipid membranes was not estab- 
lished. One should conclude from the fluorescence data that 
significant changes in the order of the lipid membrane take 
place. The evidence for phase transitions during nerve pulses 
has been discussed in more detail by Kinnunen and Virtanen 
1331 and Tasaki and coworkers I31II34I . In this context it should 
be noted that also changes in light scattering and turbidity ac- 
company the action potential that clearly cannot be related to 
membrane voltages tl4l 1351 . 



3.3 Reversible heat changes and their meaning 

The most striking thermodynamic findings in nerves during the 
action potential are reversible temperature changes and corre- 
sponding changes in the heat released during the nerve pulse. 
The first to carefully describe the heat changes was A. V. Hill 
who published a series of papers in the 1920s and 1930s. Ab- 
bott et al. 1161 showed that the heat release during the first 
phase of the action potential is nearly exactly compensated by 
a heat uptake in the second phase of the action potential. This 
effect was found in non-myelinated 1161 [T71 IT8l 1191 and in 
myelinated |16 20 1 nerves. Interestingly, Hill and collabora- 
tors found that the reversible heat release in myelinated nerves 
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originates from the complete nerve and not only from the nodes 
of Ranvier 1161 . They found it most likely that the complete 
membranes of the myelinated nerves contribute to the heat re- 
lease and that one should therefore consider an active role of 
the myelin sheet to the nervous impulse. Saltatory conduction 
that is the textbook picture for pulse propagation in myelinated 
nerves, in contrast, attributes a special role to the nodes of Ran- 
vier. Other authors reproduced the findings on reversible heat 
release, e.g. Howarth et al. 1171 . Ritchie & Keynes HI 81 or 
Tasaki and coworkers 1191 1201 l34l . It has to be acknowledged 
that these experiments are difficult and the observed tempera- 
ture changes are small (of order lOOmK). 
One important result demonstrated in Fig. 5 shows the inte- 
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Figure 5: Reversible heat change during the action po- 
tential. Left: The square of the voltage ( the energy of 
charging a capacitor) is proportional to the heat of the 
nerve pulse. The heat, however, is much larger than the 
capacitor energy. The heat during the nerve pulse returns 
to the baseline indicating that the nerve pulse is adiabatic 
( does not generate net heat after completion of the action 
potential). Data on garfish olfactory nerve adapted from 
ref. $18$. 

grated heat release during the action potential and the square 
of the voltage changes related to the free energy of the mem- 
brane capacitor 1181 . These two functions were found to be 
qualitatively nearly identical. However, the heat reversibly re- 
leased during the action potential was several times larger than 
the energy of the capacitor so that it can be excluded that the 
reversible heat release is explained by the charging of the mem- 
brane capacitor. This is the only semi-reversible element in the 
Hodgkin-Huxley model 0. Further, the heat after the whole 
pulse returns to the baseline in phase with voltage changes. 
Thus, after the nerve pulse no net heat was dissipated within ex- 
perimental error. Control experiments indicate that heat is not 



lost by thermal conduction into the environment but is rather 
reabsorbed by the nerve in the second phase of the action po- 
tential. 

The reversible heat release is a remarkable and very mean- 
ingful finding. It suggests that the physical processes underly- 
ing the nerve impulse are reversible processes. The Hodgkin- 
Huxley model, however, is based on irreversible processes, in 
particular on the exchange of potassium and sodium ions along 
ion gradients. The model does not contain any true reversible 
processes. Even if the membrane capacitor was reversibly char- 
ged, this would not result in a reversible heat change unless the 
flux of the ions was also reversible, which is not the case within 
the framework of the model. Taking the equivalent circuit pic- 
ture seriously, the flux of charges through a resistor should 
rather result in a heat release independent on the direction of the 
flux of the ions. The flux of potassium and of sodium should 
both dissipate heat. This is obviously not in agreement with the 
thermodynamic results obtained from real nerves. The finding 
of changes in lipid state and in thickness also does not find a 
satisfactory explanation within the Hodgkin-Huxley model. 

4 Propagating density pulses 

In the following we show that the thermodynamic findings de- 
scribed above find an explanation if one assumes that the action 
potential consists of a propagating density pulses. Heimburg 
and Jackson J22j showed that one could obtain stable propa- 
gating density pulses in cylindrical lipid membranes provided 
that the membrane exists in a physical state slightly above a 
melting transition. In the following we outline the underlying 
basis of this model. 

4.1 Melting transitions in biological membrane 

Many biological membranes display melting transitions slightly 
below body temperature. In Fig. 6 the melting transition of na- 
tive E.coli membranes (including all their proteins) are shown. 
One finds a pronounced lipid-melting peak slightly below body 
temperature that is affected by growth temperature of the bac- 
teria, by hydrostatic pressure and pH 1361 . Further, one finds 
several protein unfolding peaks slightly above body tempera- 
ture. It is a remarkable fact in itself that Nature chooses living 
systems to exist so close to the cooperative transitions of their 
molecules, including membranes, proteins and DNA. The un- 
derlying theme of this paper is that this is of major biological 
relevance. 

The melting transitions of such membranes display a melting 
temperature, Tm, a melting enthalpy, AH, and a melting en- 
tropy, AS, given by AS = AH/T m . Further, volume and 
area of the membrane change during the melting process. For 
the model lipid DPPC (dipalmitoyl phosphatidylcholine) that is 
the major lipid component of lung surfactant one finds: T m = 
314.2 K, AH = 35 kJ/mol, AS = 111.4 J/molK, AV/V = 



5 




temperalure |°C| 



Figure 6: Left: Schematic picture of the melting process in lipid membranes and the associated 
change in the specific heat capacity. Right: Melting profile of the membranes ofE.coli grown at 37C 
( adapted from ref. 22). The growth temperature is indicated as dashed line. The peaks below growth 
temperature belongs to the melting of lipid membranes, the peaks shaded in grey above the growth 
temperature are attributed to protein unfolding. 



0.04 and AA/A = 0.246. These values give the order of mag- 
nitude but vary between different lipid species. 

4.2 The relation between heat capacity and com- 
pressibility 

The enthalpy, specific volume and specific area changes in a 
lipid melting transition can be written as 



H(T) = H a (T) + AH(T) 
V(T) = V (T) + AV(T) 
A(T) = A {T) + AA(T) 



(4) 



H (T) is the temperature-dependent enthalpy of the pure gel 
phase and the function AH(T) is the excess enthalpy of the 
transition. Similarly, Vb(T) and A$(T) are the temperature- 
dependent specific volume and area of the gel phase. AV(T) 
and AA(T) are the excess volume and area changes associated 
with the melting transition. It has been found experimentally 
that the volume and area changes in the chain melting transition 
are proportional to the changes in enthalpy 1371 1381 . 



AV{T) = lv -AH{T) 
AA(T) = 1A -AH{T) 



(5) 



the lipid melting transition is proportional to the excess isother- 
mal volume and area compressibility: 



4(T) = 4 (T) + 4?Ac P 



4{T) 



K T,o( T ) 



V 
A 



Ac P 



(6) 



The heat capacity can easily be measured in calorimetry. The 
functions and the temperature dependent com- 

pressibilities of the pure phases that have to be taken from liter- 
ature. One can see that both volume and area compressibilities 
assume maxima at the temperature where the heat capacity is 
maximum. The adiabatic compressibilities relevant for sound 
propagation can be determined when the isothermal compress- 
ibilities are known. They assume the form 11371 



4CO 



4(T) 



T 



K 



dV\ 2 
dTj P 

2 



V-cp 
TK ' A- cp \dTj n 



(7) 



where the constants 7y = 7.8-10~ 10 m 2 /N and ^ A = 0.89 m/N 
are approximately the same for various artificial lipids and for 
biological membranes. Using the fluctuation dissipation theo- 
rem it is easy to show that excess heat capacity changes within 



where the heat capacity cp is that of the membrane plus the 
aqueous environment that transiently absorbs heat from the mem- 
brane upon compression. If the compression is very slow, cp 
will be very large and therefore in the limit of very slow com- 
pression Kg w and Kg ?s k^,. It has been found exper- 
imentally that the adiabatic compressibility obtained for peri- 
odic perturbations with a frequency ui = 5 MHz can be de- 
termined accurately if the heat capacity is assumed be the to- 
tal heat capacity of the lipid membrane alone. It is obviously 
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smaller than the isothermal compressibility. Therefore, one has 
to conclude that the adiabatic compressibility is in general fre- 
quency dependent and, thus, dispersion is present. The fre- 
quency dependence of relaxation phenomena in the lipid melt- 
ing transition has also been documented in experiments 1391 
and justified theoretically 1401 . It is also obvious from eqs. (|5j 
that the compressibility is a nonlinear function of the mem- 
brane density 1221 . If the adiabatic compressibility is known 
one can calculated the sound velocity, e.g. for the lateral sound 
velocity within the membrane plane 



K A n A 



(8) 



The lateral area density of the membrane and the enthalpy are 
related. Therefore the adiabatic compressibility is a function of 
the area density of the membrane, and it follows that the sound 
velocity is a nonlinear function of the density that, close to the 
lipid melting transition, can be expanded into a power series 
such that 



(9) 



where c is the sound velocity in the fluid phase of the mem- 
brane. Here, p and q are parameters to be determined from the 
known dependence of the sound velocity on the density. For 
unilamellar DPPC membranes slightly above the transition one 
finds experimentally that cq = 176.6 m/s (the lateral sound ve- 
locity in the fluid phase at low frequencies), p = —16.6 c^/ p A 
and q = 79.5 C 2 ,/ (p A ) 2 (for details see ref. |22)). Here, p A = 
4.035- 10 -3 g/m 2 is the lateral area density in the fluid phase of 
the membrane slightly above the melting point. Similar values 
were found for lung surfactant and native E.coli membranes. 

4.3 Propagating solitons 

We now consider the propagation of a density pulse in a cylin- 
drical membrane along the axis, x. The hydrodynamic equation 
for the propagation of such a density pulse in the presence of 
dispersion 1221 1231 is given by 



dt 2 



d_ 

dx 



ox 



(10) 



describing the changes of the lateral membrane density as a 
function of time and space. The second term is chosen ad hoc 
to mimic the frequency dependence of the sound velocity in 
a linear way using a parameter h (for details see ref. 1221 1. 
This parameter is the only one that has not yet been determined 
by experiment. We will see below that the only role of the 
parameter h is to set the linear scale of the propagating pulse. 
We have shown above that the sound velocity is a function of 
the area density, p A . Introducing eq. l|9} into eq. dlOl we obtain 



f) 2 



d_ 

dx 



(c 2 0+P (Ap A ) + q (Ap A ) 2 ) ^-Ap A 



f) 4 



(11) 
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Figure 7: Soliton profile for a soliton velocity of v=0.651 
cO calculated for h=2m4/s2. This soliton has a maximum 
amplitude of rA/rOA. Its width is approximately 10 cm. 



and after the coordinate transformation z = x—v-t (introducing 
the propagation velocity, v) we arrive at the time independent 
form describing the shape of a propagating density excitation: 



f) 2 



d_ 

dz 
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p{Ap A ) + q{Ap A Y)j-Ap' 

az 



A// 



(12) 



This equation has a localized analytical solution 1231 : 



Ap A {z) 



1 - 



v —v^ 

~ 3 J 



D" — V 
J 



■cosh ( fzjl 



(13) 

Such localized solutions are known as solitary waves or soli- 
tons. A typical soliton profile is shown in Fig. 7. The minimum 
velocity vmin allowed by eq. dl 31 is found to be 



P 
6q 



(14) 



The minimum velocity for a soliton in DPPC membranes is 
found to be v m i n = 115 m/s, which is very close to the ve- 
locity of the action potential found in myelinated nerves. The 
minimum velocity is the velocity of the soliton when its ampli- 
tude reaches the maximum value of 



A 

max 



\p\ 



(15) 
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corresponding to an overall density change of Ap^ ax /pQ = 
0.21. Soli tons with larger density change do not exist. The 
total area change when going through a melting transition is 
^Pmax/Po = 0- 246 (forDPPC). Thus, at maximum amplitude 
the soli ton forces the lipid membrane by about 85% through the 
melting transition. This will cause a transient heat release cor- 
responding to 85% of the melting enthalpy (which is on the 
order of 35kJ/mol or w 13 kT per lipid). Simultaneously, the 
thickness of the membrane will change by 85% of the thick- 
ness change in the transition from fluid to gel (7.4 Afor DPPC). 
Since the soliton is linked to changes in lipid state the fluores- 
cence anisotropy will also change. It is well known that the 
anisotropy (related to the rotational mobility) is higher in the 
gel phase then in the fluid phase. Precisely these changes have 
all been found in real nerves under the influence of the action 
potential 11411311 (see sections l3~T1 13 ,2l and l3 . 3> . The predicted 
order of magnitude of these changes matches the data found for 
such nerves. 

4.4 Electromechanical coupling 

It seems evident that the solitons described above have many 
similarities with real nerve pulses and can describe their ther- 
modynamic properties well. However, the action potential is 
known to be a propagating voltage pulse with a net voltage 
change of about lOOmV. In the following we will argue that 
this voltage change is a consequence of the change in area den- 
sity of the membrane in a manner similar to the propagation 
of a piezoelectric wave. The membranes of biological mem- 
branes contain charged lipids. Depending on cell and organelle 
the fraction of charged lipids is between 10% and 40%. Some 
membranes are especially rich in charged lipids, e.g. mitochon- 
dria. Typically, most of these charged lipids are found on the 
inner membrane, generating an electrical field. To make an es- 
timate of the size of the potential change, we therefore assume 
that the inner membrane of a nerve contains 40% charged lipids 
and the outer membrane contains only a very small fraction of 
charged lipids (average of both leaflets 20%). We ignore the 
contributions from proteins that clearly are also present. Ac- 
cording to the Gouy-Chapman theory for the potential of sur- 
faces in electrolytes, the potential of a charged surface at high 
ionic strength is given by 

* = —a . (16) 

e eK 

This is the low potential limit of the Gouy-Chapman theory 
1411 . The dielectric constant in vacuum is eo = 8.859 • 10 -12 
C 2 / Jm, and the relative permittivity e = 80 for water. Here, k 
is the Debye constant that depends on the ionic strength. For a 
monovalent salt it is given by 

/ 2e 2 

K =V^ C (17) 



where e = 1.602 ■ 10 C is the elementary charge and c is the 
concentration of the monovalent salt. For c=150mM NaCl the 
Debye constant assumes a value n = 1.26-10 9 m _1 . For a fixed 
number of charged lipids the charge density, a, is different in 
the fluid and in the gel phase of the lipids because the respec- 
tive lipid areas differ by about 24%. Therefore, one expects 
changes in the electrostatic potential of the membrane during 
a propagating density pulse. In piezoelectrics, voltage changes 
and density changes are tightly coupled. Such coupling be- 
tween lateral density and electrostatic potential is also known as 
electromechanical coupling. It is also linked to changes in ca- 
pacitance. Electromechanical coupling in membranes was first 
proposed by Petrov 1421 1431 and has been discussed by various 
authors as relevant in hair cells fl4l 1451 . Here, the potential of 
the lipid membrane is discussed. A biological membrane con- 
tains on average 50 weight percent of protein, which also carry 
charges. The total potential of the inner and outer leaflet is the 
sum of lipid and protein contributions. The contribution of the 
proteins will lead to an equilibrium resting potential of the total 
membrane that is different from that of the pure lipid mem- 
brane. However, it is most likely that only the lipids undergo 
changes in area during the pulse. The potential of the inner 
membrane at the lipid surface under the above conditions and 
the simplifying and somewhat arbitrary assumptions regarding 
lipid distribution is 

% n fi md = -114 ml/ *$ IuM = mV 

resulting in a voltage change of A^o ~ 40mV at the soliton 
peak. That is of the same order as the voltage changes in the 
action potential (which is about lOOmV). This is a very rough 
estimate since the exact charge of the lipid membrane on both 
sides of the membrane is not known and protein charges have 
not been considered. However, it seems as if the changes in 
the membrane area during the action potential are of the right 
order to account for the observed voltage changes during the ac- 
tion potential. Furthermore, membrane thickness changes dur- 
ing the action potential, thereby changing the capacitance. The 
assumption of a constant capacitance, as made by Hodgkin- 
Huxley, therefore cannot be correct (cf. eqs. and |2j . In sum- 
mary, it seems plausible that mechanical solitons can generate 
voltage changes comparable to those observed during the ac- 
tion potential. The exact values remain to be determined by 
experiment. 

5 Anesthesia 

If one assumes that the soliton model for the nerve pulse is a 
valid description of the nerve pulse containing its thermody- 
namics one immediately arrives at a quantitative explanation 
for anesthesia 1361 . Anesthesia as a tool for painless surgery 
by use of diethyl ether was first publicly demonstrated in 1846 
by William Morton from the Massachusetts General Hospi- 
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tal 1281 . This method was adopted within short time all over 
the world. Many other anesthetics had been studied in the 
following decades, including both gaseous (e.g. nitrous ox- 
ide = laughing gas) and liquid anesthetics (e.g. the alkanols 
from ethanol to decanol). A large variety of chemically dis- 
tinct molecules also cause anesthesia, e.g. barbiturates or halo- 
genated alkanes. 

5.1 The Meyer- Overton rule 

About 50 years after Morton, Meyer (53) and Overton 1261 1271 
independently found that the critical anesthetic dose of anes- 
thetics is linearly proportional to their solubility in olive oil. 
The critical anesthetic dose (or ED$o) is defined as the bulk 
concentration of anesthetic in the air (in this case equivalent 
to partial pressure) or in water at which 50% of the organisms 
studied are motionless. Overton suggested that this finding was 
related to the solubility of the molecules in the cell membrane 
whose structure was not known at the time. The Meyer-Overton 
rule covers a large range of anesthetics with membrane partition 
coefficients ranging over 5-6 orders of magnitude, from laugh- 
ing gas (N2O) and the noble gas Xenon, the liquid alcohols 
to modern anesthetics such as liducaine. The partition coeffi- 
cients of all these molecules lie within error on a straight line 
with slope 1 when plotted versus critical anesthetic dose (see 
Fig. 8, left). 

Even after more than 160 years the effect of anesthetics on 
organisms remains unexplained. A number of functions of cells 
are affected by anesthetics, including the membrane permeabil- 
ity, hemolysis, nerve function and the function of ion chan- 
nels and proteins totally unrelated to anesthesia, e.g. firefly 
luciferase. Since the most obvious effect is on consciousness 
much of the research has focused on the action of anesthet- 
ics on nerves. The Hodgkin-Huxley model is based on the 
opening and closing of ion channels, and it seems straightfor- 
ward to investigate the action of anesthetics on ion channels. In 
fact, it has been observed that some ion channel properties are 
influences by anesthetics. However, this effect is not quantita- 
tive and does not follow the Meyer-Overton rule. Some chan- 
nels are affected by some anesthetics but not by others. As 
an example, voltage gated sodium and potassium channels are 
slightly inhibited by halogenated alkanes and ethers but not by 
Xenon and nitrous oxide, although all these anesthetics follow 
the Meyer-Overton rule in causing anesthesia 1281 . It has to be 
concluded that protein pictures of anesthesia are not yet satis- 
factory. 

The Meyer-Overton rule suggests that the effect of anesthet- 
ics is independent of the chemical nature of the molecule. Since 
the noble gas Xenon lies on the same straight line as halothane 
or the liquid anesthetics, one can essentially rule out specific 
binding effects, which are the basis of the protein models (see 
also discussion). 



5.2 Melting point depression 

It is known that anesthetics have a pronounced effect on lipid 
melting transitions. Typically, with addition of anesthetics to 
the bilayers, transitions shift to lower temperatures in a lin- 
ear relation with the anesthetic concentration. Heimburg and 
Jackson 1361 have shown that this effect can be described by 
accurately by the well-known phenomenon of freezing point 
depression. If one assumes that anesthetics molecules are read- 
ily soluble in fluid lipid membranes and insoluble in the gel 
membrane, one arrives at following law for the freezing point 
depression 

ATm = ~^H XA ' (19) 
where xa is the molar fraction of anesthetics in the fluid lipid 
membrane, T m is the melting point of the lipid membrane and 
H is the melting enthalpy. The derivation of this equation can 
be found in any physical chemistry textbook. The membrane 
concentration of anesthetics at the critical dose, xa, is related 
to the partition coefficient via 

x A = P ■ {ED 50 ) ■ Vi (20) 

where P is the partition coefficient between membrane and 
water, -E-D50 is the critical anesthetic concentration, and Vi is 
the molar volume of the lipids (about 0.75 1/mol). The two 
equations above describe the behavior of many anesthetics. D. 
Kharakoz |46l has collected data for various anesthetics, some 
of which are displayed in Fig. 8 (right). Shown is the concen- 
tration dependence of the melting point as a function of the crit- 
ical anesthetic dose for tadpoles. The melting point depression 
for all anesthetics (shown here are alkanols) lie on a straight 
line when plotted versus the critical anesthetic dose. The slope 
of the curve indicates that the shift of the transition temperature 
at critical anesthetic dose is AT m = —0.6 K for all anesthetics 
that follow the Meyer-Overton rule, independent of the chemi- 
cal nature of the drug |36 46 1. The Meyer-Overton rule there- 
fore can be reformulated as: The anesthetic potency of anes- 
thetics is proportional to their ability to lower the melting point 
of lipid membranes. It is clear that within the soliton model 
for nerve pulses the melting points play an essential role. The 
assumption in the following is that the lipid melting point plays 
an important role in the control of biological membranes. 

5.3 Pressure reversal 

If one assumes that the lipid melting point is important for bi- 
ological function and that the effect of anesthetics is related to 
their effect on melting points, it is interesting to compare this 
to other physical properties that also influence melting points, 
most notably the influence of pressure. It has long been known 
that pressure influences the melting points of membranes by 
shifting them to higher temperatures. The pressure dependence 
of such transitions is described by 1381 

AT m = 7v ApT m (21) 
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Figure 8: Left: The Meyer-Overton rule for volatile anesthetics showing the linear dependence of 
the oil/gas partition coefficient and the critical anesthetic dose for man. The solid line represents a 
straight line with slope 1. Data adapted from ref. 27. Right: Lowering of the melting transition for 
a series of alkanols as a function of the critical anesthetic dose for tadpoles. The solid line displays 
a slope of 1. Adapted from ref. 46. 



where -fv = 7.8 • 10 _1 m 2 /N is a constant that is roughly the 
same for all lipids, lipid mixtures and biological membranes 
1371 1381 . This equation indicates that a lipid membrane with a 
T m = 314 K (dipalmitoyl phosphatidylcholine) shifts its transi- 
tion by 1 K to higher temperatures upon application of 40.8 bar 
hydrostatic pressure. This indicates that a pressure of 24.5 bar 
should be sufficient to reverse the effect of anesthesia (which 
corresponds to a shift by 0.6 K to lower temperatures). Pressure 
reversal of anesthesia has indeed been found, first by Johnson 
et al. 1471 . If tadpoles are anesthetized at 3 times the criti- 
cal anesthetic dose of ethanol, they wake up upon application 
of 150 bars of hydrostatic pressure. The pressure reversal of 
anesthesia is well documented in the literature. 

5.4 Free energy of the membrane 

The free energy difference between gel and fluid phase is the 
free energy that must be provided to shift the lipid membrane 
through its phase transition. It is given by 



AG = AH — TAS = AH 



T m -T 

T 



(22) 



making use of the identity AS = AH/T m . This equation indi- 
cates that the free energy difference between the two phases is 
linearly dependent on the difference of the experimental tem- 
perature, T, and the melting temperature, T m . Now, we have 
shown in the previous section that T m is influenced by both 
anesthetics concentration and by pressure. The melting tem- 
perature T m is changed by anesthetics and pressure in the fol- 
lowing manner 



T — T 



m,0 



771 

Ag XA + IvApTrnfi (23) 



where T m o is the transition temperature at atmospheric pres- 
sure and in the absence of anesthetics. We finally obtain 



AG(x A , Ap) w AH 



T, 



rn.O 



- T 



RT . T 
- ~^ X A + lv Ap- 



(24) 

If the melting transition of the lipid membrane is to play a rele- 
vant role for biological function, it follows that biological func- 
tion should be the same when AG is the same. Therefore, the 
condition for pressure reversal of anesthesia is or 



Ap 

APED 50 



1 RT m ,o 
^^H~ XA 
1 RT, 

lv 



ih^Lp^ED^. (25) 



The numbers obtained from this equation are of an order very 
similar to that obtained in experiments. Data from octanol and 
DPPC membranes as well as the equations above suggest that 
a pressure of 24.5 bar reverses anesthesia 1361 . The data from 
Johnson on tadpoles in an ethanol solution corresponding to 
three times the anesthetic dose was reversed by 150 bars of 
pressure 1471 . Our calculation yields 73.5 bars, assuming a 
membrane partition coefficient for ethanol of 0.6 (which is sub- 
ject to an error of the order of a factor 2). 

6 Discussion 

We have suggested here that the Hodgkin-Huxley model (6) for 
the action potential does not provide a satisfactory description 
of nervous impulse because it does not include the mechanical 
and optical changes associated with the action potential. Fur- 
ther, it is clearly inconsistent with the thermal response. The 
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initial heat release and subsequent re-absorption studied by a 
number of authors 1161 [T71 IT8l |T9l 1201 points rather to a re- 
versible physical phenomenon that conserves entropy. In con- 
trast, the Hodgkin-Huxley model is based on the flux of cur- 
rents through resistors that should heat the membrane indepen- 
dent of the nature of the current and its direction. Therefore, 
we have proposed an alternate model based on the known me- 
chanical and thermal features of artificial and biological mem- 
branes. It was shown that under physiological conditions sta- 
ble mechanical solitons could propagate and display reversible 
heat release, changes in membrane thickness, changes in mem- 
brane order and reversible membrane potential changes. All 
these changes have been observed in experiments. In particu- 
lar the reversible heat release and the overall conservation of 
entropy is a feature typical of sound propagation. It should be 
noted that we use the term sound propagation in a general sense 
that includes all changes of the thermodynamic variables that 
accompany a mechanical compression according to Maxwells 
relations. In such a description, the simultaneous occurrence of 
density changes, voltage changes, and heat release is a surprise 
but rather a necessary consequence of thermodynamics. The 
Hodgkin-Huxley model seems to be in agreement with fluxes 
through ion channel proteins. However, the currents through 
such channels fall short of presenting an explanation in the 
sense of a physical theory based on first principles. The con- 
ductances of the channels contain many parameters that can- 
not be justified theoretically. Therefore, their seemingly simple 
description relies on objects that contain all the unexplained 
features in the form of parameters. For this reason Hodgkin 
and Huxley originally recommended treating their model with 
care. They state in their seminal paper from 1952 |6|: "The 
agreement must not be taken as evidence that our equations 
are anything more than an empirical description of the time- 
course of the changes in permeability to sodium and potas- 
sium. An equally satisfactory description of the voltage clamp 
data could no doubt have been achieved with equations of very 
different form, which would probably have been equally suc- 
cessful in predicting the electrical behavior of the membrane. 
. . . the success of the equations is no evidence in favour of the 
mechanism of permeability change that we tentatively had in 
mind when formulating them." In this paper we have, in fact, 
shown that many changes can be explained by totally differ- 
ent physical mechanisms that result in similar equations for the 
pulse propagation. Hodgkin was clearly aware of the problems 
generated by the finding of a reversible heat releases during the 
action potential. He wrote in his textbook 'The conduction of 
the nervous impulse' 1211 : "In thinking about the physical ba- 
sis of the action potential perhaps the most important thing to 
do at the present moment is to consider whether there are any 
unexplained observations which have been neglected in an at- 
tempt to make the experiments fit into a tidy pattern. . . .perhaps 
the most puzzling observation is one made by A.V. Hill and his 
collaborators Abbott and Howarth ( I958). Xl6\l . . . Hill and his 
colleagues found that it (the heat release) was diphasic and 



that an initial phase of heat liberation was followed by one of 
heat absorption, a net cooling on open-circuit was totally un- 
expected and has so far received no satisfactory explanation." 
Howarth et al. 1171 concluded from their finding of heat re- 
lease and subsequent heat uptake: "It seems probable that the 
greater part of the initial heat results from changes in the en- 
tropy of the nerve membane when it is depolarized and repolar- 
ized." Reversible entropy changes, however, are not a feature 
of textbook pictures of nerve pulses. Here, we have followed 
Hodgkins suggestion and searched for ways to explain the re- 
versible heat. Slightly below physiological temperatures, there 
exist chain-melting transitions of the membrane. It is interest- 
ing to note that these transitions occur at much lower tempera- 
tures in the absence of the proteins. For instance, the melting 
point of E.coli lipid extracts is about 20K lower than that of 
the native membrane in the presence of all their lipids. There- 
fore, the presence of proteins seems to play an essential role in 
fine-tuning the thermodynamics of biological membranes. Be- 
sides their role as catalysts, proteins also possess chemical po- 
tentials that are thermodynamics variables. They contribute to 
the behavior of membranes in a manner similar to temperature, 
pressure, pH and other variables. The presence of cooperative 
lipid transitions forms the basis for the possibility of density 
pulses that propagate along the nerve axon. One short-coming 
of our model is that it does not yet include a frictional term 
even though one may expect that, due to the flux of lipids and 
changes in diameter of the nerve, a proper hydrodynamic treat- 
ment should yield in a dampening of the pulse. This problem 
remains unanswered in the context of our model, mainly due to 
the lack of detailed data on the dimensional changes in nerves. 
However, experiments show that such density pulses propagate 
l9l llOI in real nerves, and the near-complete reversal of the heat 
ll6M17lfT8l suggests that friction is small. Within the soliton 
model proteins do not play a role as channels or as active com- 
ponents. Rather, they tune the thermodynamics of the mem- 
brane. An important question is how such a mechanical soli- 
ton can be generated in a membrane. Since the soliton pushes 
the membrane through its chain melting transition, everything 
that moves membranes through transitions should be able to 
generate a pulse. All physical changes that push the transition 
away from physiological conditions should inhibit pulses. As 
an example, local cooling of a nerve has been shown to induce 
nerve firing, whereas temperature increase inhibits pulse con- 
duction 1311 . Due to the electromechanical coupling described 
in section 14741 changes in trans-membrane voltage are also po- 
tentially able to generate pulses. Further, a local decrease of 
pH, increase in pressure or increase in calcium concentration all 
have the potential to trigger pulses because all of these changes 
increase the phase transitions of biomembranes. Most interest- 
ingly, anesthetics inhibit pulse generation due to their property 
of lowering phase transitions. Since ion channels do not play 
an active role in our description of nerve pulses, it is obvious 
that the action of anesthetics requires a different explanation 
than their action on ion channels. The famous 100-year old 
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Meyer-Overton correlation 1251 l26l 127 28 1 states that the ac- 
tion of anesthetics is, within error, strictly proportional to their 
solubility in lipid membranes. This law is valid over 6 orders 
of magnitude in the membrane/air and membrane/water parti- 
tion coefficient. This law remains an elegant and valid means 
to determine the effectiveness of an anesthetic 1281 . It basically 
excludes the notion that the action of anesthetics can be linked 
to specific binding of the drug to a receptor. The argument is 
simple: The binding of two molecules is described by the free 
energy, which is a function of state. If the action of anesthet- 
ics is exactly proportional to the concentration of drugs in the 
membrane independent on chemical nature of the drug as fol- 
lows from the Meyer-Overton correlation, the binding constant 
of all anesthetics to receptors must be identical, including that 
of the noble gas Xenon. Since noble gases cannot bind specifi- 
cally, the same must be concluded for all other anesthetics that 
follow the Meyer-Overton rule. The experimental finding is 
that halogenated alkanols act very differently on ion channels 
than Xenon or nitrous oxide. Thus, protein models are clearly 
not consistent with the well-documented Meyer-Overton cor- 
relation. Although protein models are currently quite popular, 
they cannot fulfill the basic thermodynamic requirements for 
anesthetics that follow the Meyer-Overton correlation. Due to 
the above argument it is unlikely that the action on ion channels 
is related to anesthesia. Here, we have outlined the thermody- 
namic theory of how anesthetics influence the phase behavior 
of lipid membranes via a well-known unspecific phenomenon 
known as freezing-point depression. It states the lowering of 
the melting point is proportional to the membrane concentra- 
tion of the anesthetic drug. Thereby, we attribute a physical 
meaning to the Meyer-Overton rule that was not provided by 
Overton himself. By this mechanism anesthetics alter the fea- 
tures of propagating solitons in a quantitative manner. More 
specifically, they alter the amount of free energy that has to 
be provided to generate a pulse. We found that it is linearly de- 
pendent on the distance between physiological temperature and 
the transition in the nerve membrane. This approach admits the 
possibility of finding strict thermodynamics relations between 
various thermodynamics variables, including the pressure re- 
versal of anesthesia that can be calculated in quantitative terms 
1471 . It seems unlikely that all of these quantitative correlations 
can be found experimentally without thermodynamics being an 
essential player in the description of both the action potential 
and the action of anesthesia. Indeed, simple thermodynamics 
seems to contain a complete description of such phenomena. 
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